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We study the q-dependent susceptibility x(q) of a Z-invariant ferromagnetic 
Ising model on a Penrose tiling, as first introduced by Korepin using de 
Bruijn's pentagrid for the rapidity lines. The pair-correlation function for 
this model can be calculated exactly using the quadratic difference equations 
from our previous papers. Its Fourier transform xil) is studied using a novel 
way to calculate the joint probability for the pentagrid neighborhoods of the 
two spins, reducing this calculation to linear programming. Since the lattice 
is quasiperiodic, we find that x(q) is aperiodic and has everywhere dense 
peaks, which are not all visible at very low or high temperatures. More and 
more peaks become visible as the correlation length increases — that is, as 
the temperature approaches the critical temperature. 

KEY WORDS: Ising model; quasiperiodicity; Fibonacci sequence; pentagrid; 
Penrose tiles; Z-invariance; correlation functions; q-dependent susceptibility. 



1. Introduction 

In an experiment® done in 1984, Shechtman and his coworkers found fivefold 
symmetry in the diffraction patterns of some rapidly cooled alloys. As such a 
symmetry is incompatible with lattice periodicity, it was concluded that the 
crystalline structures of these alloys, if any, must necessarily be quasiperiodic. 
This theoretical explanation came forward almost immediately, as Penrose, 
de Bruijn, and Mackay'l^'^'^'^'SlO had already studied tilings that have fivefold 
symmetry, well before this experimental discovery. Quasiperiodic tilings are 
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types of almost periodic structures that permit sharp peaks in the diffraction 
patterns, but have normally forbidden symmetries. ''^'^'^'^'^'^ 

Already in 1986, Korepin'El' introduced a Z-invariant eight- vertex model 
on Penrose tiles. The Z-invariant inhomogeneous models are completely in- 
tegrableGS even on irregular lattices and their critical exponents are known 
to be the same as those of homogeneous systems on regular lattices. Thus, the 
critical behaviors of these quasiperiodic Z-invariant mo dels C31l511I5i n^[i9l E5]|2iJllIli 
have to be the samejfl independent of the lattice structurejll 

In such models, the order parameter is the same''^'^ for all sites and it 
vanishes towards the critical point. Therefore, the Fourier transform of the 
one-point function of a Z-invariant Ising model is the product of this order 
parameter and the lattice sum ^ e''*''". Experiments that probe the resulting 
"magnetic" Bragg peaks are restricted to the low-temperature phase and the 
corresponding theory is essentially the zero-temperature theory, well-studied 
in the literature. ''^'^'^ The aforementioned Bragg peaks will broaden, if we 
allow the underlying quasicrystalline lattice to become distorted by lattice 
vibrations . I^l^ However, we shall not consider this possibility in this paper, 
as we assume the underlying lattice to be a perfect and rigid Penrose tiling, 
restricting our attention solely to the ordering of the spins at the lattice sites 
under thermal fluctuations. 

Contrary to the above theory for the Bragg scattering, the situation is far 
more complicated for scattering experiments that probe the pair-correlation 
function (o"rCrr') via the wavevector-dependent susceptibility x(q). This last 
quantity is defined as 

fcBTx(q) = X(q) = jim ^ e''^'^^'-) [(arCXrO - {a,)M], (1.1) 

r r' 

where C is the number of lattice sites, r and r' run through all these sites, 
and q = {qx, Qy). It is the Fourier transform of the connected pair correlation 
function, defined by what is inside the square brackets. Furthermore, x(q) is 

Universality of the critical exponents of ferromagnetic Ising models on quasiperiodic 
lattices has been confirmed for non-Z-invariant cases also using real-space renormaliza- 
tion group techniques, ''23 Monte Carlo simulations, ''^'^'^UlSllSZllSSlllg) gej.jgg expansion 
methods, <3Sl[3ll|33 ^nd the study of Yang-Lee zeros-^^SMm 

Unlike the q-dependent susceptibility, thermodynamic quantities like the free energy, 
the specific heat, and the bulk susceptibility do not probe the lattice structure, although 
subtle lattice effects do show up in corrections to scaling. 'SSJ 
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the reduced q-dependent susceptibility, taking out a trivial factor involving 
the absolute temperature T. 

In the well-known lattice-gas language it becomes proportional to the 
structure function, the Fourier transform of the density-density correlation 
function, which is also measurable in diffraction experiments, revealing the 
symmetry of the lattice. Thus, the ^'-invariant models can indeed be 

used to show the difference between quasiperiodic and regular lattices and 
we expect it to provide a diffuse scattering pattern both above and below 
the critical temperature Tc, with more structure closer to Tc. 

There is good reason to pick a Z-invariant Ising model for our present 
study of It is taken from the foremost class of models with short-range 

interactions allowing exact computations. Moreover, comprehensive exten- 
sive studies on spin-spin correlation functions in nontrivial models with short- 
range interactions have been done only for Ising models. '^^^ '^^ 43 ,44 ,45 46 ,47, iUHHli 
Cited here is just a fraction of the literature. The accumulative knowledge 
of these studies has made the calculations in Z-invariant quasi-periodic Ising 
models possible. Correlations in other nontrivial cases are still mostly inac- 
cessible to exact methods of evaluation. 

For instance, in 1988, Tracy'''^^^^' introduced the layered Fibonacci Ising 
model, which is not Z-invariant. The row correlation functions in the layered 
Ising model are known to be block Toeplitz determinants,''^ of which we still 
have no idea how to evaluate them exactly in general, except in a few simpler 
cases. ''^'^ Tracy has shown that the critical exponent of the specific heat 
remains unchanged in such quasiperiodic layered models. 

In our previous papers, ''^'^ we have studied Fibonacci Ising models 
whose spins are on regular lattices, but whose nearest-neighbor interactions 
are quasi-periodic. They are special cases of inhomogeneous Ising models 
whose Hamiltonians are given by 



with /3 = 1/k-BT. When the system is periodic, the pair correlations are 
translationally invariant. Thus one of the sums in (11.11) can be carried out. 
For T ^ Tn, the connected correlations decay exponentially as functions 
of distanced Therefore, only a finite number of short-distance correlation 
functions are needed for the calculation of the iii Eq- (H-IP - 

^ For T — Tc, the correlations decay algebraically, so that we need the Epstein-Ewald 
summation formula to take into account the long-distance behavior. 



pn=J2 (Km, 
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For the Ising models considered previously, ''^'^ the couplings between 
specific nearest-neighbor spin pairs form Fibonacci sequences {Sn} defined 
recursively''^ by 

Sn+l = SnSn~l, Sq = B, Si = A, (1-3) 

so that 5*2 = AB, S3 = ABA, S4 = ABAAB and so on. Since this sequence 
is quasiperiodic, as arbitrary long subsequences are repeated infinitely often, 
the model is also aperiodic. Consequently, the correlations are no longer 
translationally invariant. However, the averages of the correlations for two 
spins at fixed distance can be evaluated by using a theorem of Tracy. 'S3' 

In our previous works, ''^'^ we have found that the various q-dependent 
susceptibilities of our Fibonacci Ising models are always periodic. They 
can have multiple incommensurate everywhere- dense peaks in each unit cell 
only if the aperiodic oscillations in the average correlation functions are not 
negligibly small. This is true in the mixed case when the interactions are 
aperiodic sequences of ferromagnetic and antiferromagnetic couplings. The 
number of visible peaks in increases as the correlation length increases. 
In contrast, in the periodic Fibonacci Ising lattices with mixed bonds, the 
visible peaks in the q-dependent susceptibility are at commensurate positions 
and their number has a finite maximum.''^ 

The ferromagnetic aperiodic Fibonacci Ising lattice, on the other hand, 
behaves almost like the regular Ising model — one peak per unit cell, located at 
the commensurate position — because the aperiodic oscillations in its average 
correlation functions are negligibly small at all temperatures.''^ 

The q-dependent susceptibility is found to have almost the same 

behaviors for both T < and T > T^.^ This is in agreement with the result 
of Peter Stephens, who showed that the randomized (disordered) icosahedral 
system''^ gives almost the same diffraction pattern as a quasicrystal — which 
is in the solid (ordered) phase. 

Even though both ferromagnetic and antiferromagnetic edge interactions 
are present in the mixed case, the mixed systems in Refs. and EH are 
not frustrated. In fact, their partition functions are equal to the partition 
functions of the ferromagnetic models, from which they differ by gauge trans- 
formations of signs. Also, the x(q) of a fully-frustrated model may not show 
incommensurate peaks. ''^ 

In this paper, we turn our attention to systems with a quasiperiodic lattice 
structure. H^JIMIISD More specifically we study a Z-invariant Ising model whose 
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spins are on vertices of a Penrose fat-and-skinny rhombus tiling with the 
pentagrid as its rapidity hnes.''^'!^'!^'!^'!^ 

1.1. Outline 

This paper is organized as follows. In Section [2], we introduce a Z-invariant 
Ising model on Penrose tiles constructed from a pentagrid,'^' with spins on 
either the odd or even sublattice. We show how the pair-correlation functions 
can be evaluated in Section [21 Section H] is rather lengthy. In it, a new 
method of counting all the spin sites and evaluating the joint probabilities 
of the occurrence of two neighborhoods of two spins is given. In subsection 
15.11 we describe in detail the calculation for the q-dependent susceptibility 
for the odd lattice. The results are given in subsection 15.21 In subsection 15.31 
we present a mapping between the odd and the even sublattices. Finally, in 
Section [6] we present our conclusions. 

2. Pentagrid and Penrose Tiles 

In two ingenious papers by the famous Dutch mathematician N.G. de Bruijn, 
he relates the non-periodic Penrose tilings in a plane to a pentagrid,® which 
is a superposition of five grids. Each grid consists of parallel lines with equal 
spacings between the lines; the grids may be obtained from one another by 
rotations of angles which are multiples of 27r/5. This is shown in Fig. [1] Here 
some grid lines of the pentagrid are shown; the arrows on the lines shown in 
Fig. [1] should be ignored for the moment, as they will define the direction of 
the "rapidities" which we define later. 

To describe the pentagrid in mathematical formulas,® let 

^ = e^Ws, C + C"' = 2 cos(27r/5) = = i( - 1), (2.1) 

in which p is the golden ratio. Then, choose 7o, 7i, 72, 73, 74 to be five real 
numbers, satisfying 

7o + 7i + 72 + 73 + 74 = 0. (2.2) 
Now the jth grid in the pentagrid consists of lines given by 

Gj = {z E C|Re«--') + 7j- = kj, kj G Z}, j = 0, ■ ■ ■ , 4. (2.3) 

The pentagrid is called regular, if there is no point in the complex plane C 
belonging to more than two of the five grids. This also means, every vertex 
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Fig. 1. The pentagrid is a superposition of five grids, each of which consists of parallel 
equidistanced lines. These grid lines are the five different kinds of rapidity lines in a Z- 
invariant Ising model. Arrows are put on the grid lines, and all of them pointing to the 
upper half plane. Through a center in one of the meshes (site for a spin cr), dashed lines 
are drawn to split the pentagrid into ten regions I to X. Correlations between o and a' 
are different when a' sits in different regions. Five little arrows perpendicular to the five 
dashed lines indicate the directions in which the corresponding integers fcj, j = 0, • • • , 4, 
increase. 
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of the regular pentagrid is an intersection of no more than two hnes. Each 
vertex is surrounded by four meshes (which are often called faces in physics) . 

Now to every point z in the complex plane C, de Bruijn associates an 
integer vector K{z) = {Ko{z),- ■■ ,K4{z)) whose five elements are integers 
given by 

K^{z)=\Re{zC')+^,], (2.4) 

in which [x] denotes the "roof of a;", which is the smallest integer ^ x. 
It is easily seen from 02.41) and (12.31) that whenever z moves across a line 
of the jth grid, Kj{z) changes by 1. All points in the same mesh (face) 
have the same integer vector and the integer vectors of different meshes are 
different. From Fig. [1], we may already see that some fraction of the meshes 
(or faces) becomes infinitesimally small in size as the number of lines in each 
grid becomes infinite. 

Since every vertex of the regular pentagrid is surrounded by four meshes 
(faces), by assigning to each of their four corresponding integer vectors K{z) 
= {Ko{z), ■ ■ ■ , K^lz)) a complex number 

4 

f{z) = J2KA^)C, (2.5) 

j=0 

these four meshes are now mapped to the vertices of a rhombus. More 
specifically, to the intersection of two grid lines kr and kg, (r 7^ s), one 
assigns a rhombus in C whose vertices are the four complex numbers f{z), 
f{z) + f{z) + and f{z) + C + C'^ assigned to the four surrounding 
meshes. Clearly, there are two different kinds of rhombuses: the thick one 
having angles 72° and 108° for r = s ± 1, and the thin one having angles 36° 
and 144° for r = s ± 2. In both rhombuses, all sides have length 1. They are 
shown in Fig. [2l 

In these two papers,® de Bruijn also showed that, even though there are 
many different choices of 7^ in (12.21) . many of the resulting pentagrids are 
shift-equivalent, that is, they can be obtained from each other by a parallel 
shift @ 

We now assign to each grid line in the grid j, j = 0, ■ ■ ■ , 4, of the pentagrid 
a rapidity uj pointing into the upper half plane, as is shown in Figs. [1] and 

^ He also proved that if ^ = 7jC^'' times a power of C is not purely imaginary (modulo 
the principal ideal of 1 — C given by all complex numbers of the form ^jC"' with integers 
TiQ, • • • ,n4 satisfying ^ rij — 0), then the corresponding pentagrid is regular.® 
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(a) Thick rhombus (b) Thin rhombus 



Fig. 2. Spins interacts along the diagonals. All sides have length 1. The long diagonal 
of the thick rhombus has length p and the short diagonal of the thin rhombus has length 
p~^, with p the golden ratio, given in (j2.1|) . The little arrows perpendicular to the grid 
(rapidity) lines indicate the directions in which the integers kj, j = 0, - ■ ■ ,4, increase. 

[3l In Fig. [HI it is also shown how the five grid fines kj = 0, for j = 0, ■ ■ ■ , 4, 
are sfiifted to make tfie grid regular, again following de Bruijn.® 

Tfie usual Z-invariant Ising model''^SEli jg formed by putting spins inside 
tfie mesfies, but fiere, fiowever, tfie Ising spins are on tfie vertices of tfie 
rfiombuses. Tfiere is an one-to-one mapping given by (12. 5p relating a mesfi 
in tfie pentagrid to a vertex of tfie Penrose tiling. After tfie mapping, tfie grid 
lines in tfie pentagrid — wfiicfi are also tfie rapidity lines — become "Conway 
worms" (no longer straigfit) in tfie Penrose tiling.''^ Since tfie rapidity lines 
are used to define commuting transfer matrices, tfiey do not fiave to be 
straigfit lines. 

Tfie model so defined is a special case of tfie infiomogeneous Z-invariant 
eigfit- vertex model proposed by Korepin. USJlisiCSMIB More specifically, tfie 
four-spin couplings are identically zero, and tfie eigfit-vertex model decom- 
poses into two independent Ising models. Tfie interactions of tfie Ising spins 
are along tfie diagonals of tfie rfiombuses. Tfie odd and even sublattices are 
tfierefore decoupled. ''^'^ 

As in tfie earlier works, ''^'^ tfie coupling K{ui, vj) between two spins is 
represented by a line connecting tfiese two spins, witfi tfie arrows of tfie two 
rapidity lines Ui and Vj on tfie same side of tfiis line, as sfiown in Fig. 111(a), 
wfiile tfie line representing tfie coupling K{ui, Vj) fias tfie arrows of tfie two 
rapidity lines Ui and Vj on opposite sides, as sfiown in Fig. Il](b). Tfie edge 
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Fig. 3. The five different kinds of rapidity lines uq, ■ ■ ■ ,1*4 pointing into tlie upper lialf 
plane. The grid lines kj = 0, for j = 0, • • • ,4, are shifted from the dashed lines to make 
the pentagrid regular. 

interactions are parametrized by 

sinh (2K(ui,Vj)) = ksc{ui — Vj, k') = cs(A + Vj — Ui, k'), 

sinh (2K{ui,Vj)^ = cs{ui — Vj,k') = k sc(^X + Vj — Ui,k'y (2.6) 

Here A = K{k') is the elliptic integral of the first kind, and k and k' = yX^^W 
are the elliptic moduli; these are convenient temperature variables, assumed 
to be the same for all sites. 



'J 



(a) 



(b) 



Fig. 4. (a) The horizontal coupling K{ui,Vj); and (b) the vertical coupling K{ui,Vj). 
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From Fig. O we can see that the lengths of the four diagonals of the 
two rhombuses are different. The interactions between the spins are chosen 
to depend on the interparticle spacings only, but not on the orientations. 
Consequently, we must have 

Uq — Ui = U2 — Us = U4 — Uq = X + Ui — U2 = X + U3 — U4. (2-7) 

From this, we find 

4A 3A 2A A ,^ 

U4-Ui = —, U2-Ui = —, Uo-Ui = —, U3-Ui = -. (2.8) 

If we let 

Sj = ksc{jX/5,k'), (2.9) 

then for the thick rhombus in Fig. [2](a), we assign S2 to the longer diagonal 
and S3 to the shorter diagonal, while for the thin rhombus in Fig. [2](b), S4 
to the shorter diagonal and si to the longer one. Thus, to the four types 
of diagonals are assigned four kinds of couplings according to their lengths, 
with a stronger coupling for a shorter interparticle distance. 

The two Ising sublattices on the Penrose tiling are indicated in Fig. O 
The edges in the even sublattice are omitted. There are eight types of vertices 
S, K, Q, D, J, S3, S4, S5 in the Penrose tiling, which are shown in Fig. 7 of 
Ref. O The coordination numbers of spins in the Penrose Ising model are 3 
for types Q and D; 4 for K; 5 for S, J and S5; 6 for S4; 7 for S5. 



3. Correlations 



Two spins (Tr and ar' at two different vertices of the Ising lattice just defined 
have different integer vectors K = {Kq, ■ ■ ■ , K4) and K' = {K'q, ■ ■ ■ , K'^). 
Since there is a one-to-one mapping between the vertices of the Penrose tiles 
and the meshes of the pentagrid, the integer vectors can be used to denote 
the positions of the spins: r ^ K. 

Since each grid consists of parallel lines with equal spacings, the absolute 
value of the difference ij = Kj — Kj is actually the number of the jth kind of 
rapidity lines sandwiched between these two spins. The correlation functions 
were 

shownlHEi to be 

(o-^cr^,) = (aa')[^o,...,£4] 

\io\ \i2\ lisl K^l 

~ fi'(^0' • • • ' ^0' ^1' • • • ' ^1' ^2' • • • ' ^2' ^3' • • • ' ^3' ^4' • • • 5 ^4)5 (3-1) 
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Fig. 5. Two Ising lattices on Penrose tiles. The spins on the odd sublattice (the first 
Ising model) are joined by bonds, but the spins on the even sublattice (the second Ising 
model on the dual lattice) are denoted by circles without their bonds. The two lattices are 
independent. The coordination number, which is the number of adjacent spins connected 
to a given spin in each sublattice, takes values 3, 4, 5, 6, or 7. Incorporating the two models, 
a third Ising model can be defined with Ising spins on both the odd and even sublattices, 
i.e. the decoupling limit of Baxter's eight-vertex model on the Penrose lattice . GU P^Iisl lISlllH 
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where u'j = uj for rapidity lines of type j with arrows pointing to the same 
side of the hne joining the two spins, and u'j = uj ± X for rapidities with 
arrows pointing to opposite sides of the hne. It is as if the rapidities of hnes 
pointing to the opposite side need to be flipped by adding ±A, i.e. adding ±n 
to the angle variable Ujir/X.^ The functions g have both the "permutation 
sjTumetry" (which means that they are invariant under all permutations of 
the rapidities) and the "difference property" (which implies a translation 
invariance under shifting all the rapidities by the same amount).'^ 

We next examine in more detail when we have to choose u'j = Uj or 
u'j = Uj±X in (13.11) . Through a point in the mesh where spin a sits, as shown 
in Fig. [H we draw five dashed lines parallel to each of the five grids. As the 
choice of the point in the mesh is rather arbitrary, the dashed lines should 
really have been drawn with a finite thickness, i.e. open strips without their 
boundaries between two consecutive grid lines of the pentagrid. Five arrows 
are also drawn perpendicular to the dashed lines to indicate the directions 
in which the integers kj increase. If the other spin a' is in a mesh crossed 
by the jth dashed line, i.e. a and a' lie within the same strip, then the two 
spins have the same Kj {ij = 0) and their pair correlation function does not 
depend on the value of Uj. When a' moves away from this dashed line in the 
direction of the arrow, we have ij ^0, whereas ij ^ if a' moves away in 
the direction opposite to the arrow. These dashed lines (or more precisely 
strips) divide the entire plane into ten regions, and we numbered them from 
I to X. From Fig. [1], we can see that [io, ■ ■ ■ , ^4) have the same signs inside 
each region. 

If a' is in regions III or VIII, we find from Fig. [1] that the arrows of all 
the rapidity lines point to the same side of the line joining the two spins. 
Thus u'j = Uj for all five j-values. When a' is in regions II or VII, then the 
rapidity lines with M4 and the other rapidity lines are pointing to opposite 
sides of the hne joining the spins, so that u'^^ = U4 — X and u'j = Uj for j 7^ 4. 
For regions IV and IX the Ui rapidity lines point in the other direction with 
respect to the other rapidity lines, implying u[ = Ui + X. If a' is in regions 
I and VI, the arrows of the U2 and U4 rapidity lines are on the opposite side 
and, therefore, = U4 — A and = M2 — A. Similarly, for regimes V and X, 
u[ = Ui + X and = M3 + A. For all other j-values, u'j = Uj. In summary, 
our choices for u'j are listed in Tabled with the Uj's given in (12.81) . where we 
may set Ui = without loss of generality in view of the difference property 
of the pair correlation function. 

We may even shift the five rapidity values Mq, . . . , U4 in (13.11) by the same 
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Table I. Rapidities 



Regions 


U4 


U2 


^0 


^3 


u[ 


I and VI 


M4 — A 


U2- X 


Uq 


U3 


Ui 


II and VII 


M4 — A 


U2 


Uo 


U3 


Ui 


III and VIII 


U4 


U2 


Uo 


U3 


Ml 


IV and IX 


U4 


U2 


Uq 


U3 


ui + X 


V and X 


U4 


U2 


Uq 


U3 + X 


ui + X 



amount, depending on the choice of region, such that min^w^ = 0. We can 
then also use the permutation property''^ of the pair-correlation function g, 
given in fl3.ip . to rearrange the five resulting rapidity values u'j in decreasing 
order as |A,|A,|A,|A,0. Therefore, it suffices to calculate the quantity 

g[m4,, 7713, 1712, mi, mo] = 



m,4 



1713 



m2 



mi 



mo 



( 4A 4A 'SA 



3A 2A 



2A A 



A 



■' 5 ' 5''"' 5'°' 



.0 



(3.2) 



where the m^-'s are nonnegative integers depending on the ^j's and the choice 
of region. To determine this dependence, we can first use (12. 8p . from which 
we find M3 + A > Ml + A > M4 > ^2 > Mo > > Ml > M4 — A > M2 — A. 
Comparing (13.11) . Tableland (12.81) we can then express all pair correlations 
in the form (13.21) . We list the results for the ten different regions in Table HTl 
This completes, more or less, the calculation of the pair correlation function, 
as we can refer to our previous papers''^'^ for further details on how to 
evaluate the g\mii,m3,m2,m\,m^ defined in (13. 2p . 

Since the pentagrid is invariant under rotations by angles that are integer 
multiples of 27r/5, the grids may be relabeled m m + j , (mod 5). Then 
the differences of the five integer vectors of the two spins are also relabeled 
^m ^m+j- This shows that the pair correlation function must have the 
cyclic property 



(o-O-')[£0A,;?2/3/4] = ('^'^')[^,/i+l/i+2/i+3/i+4]5 (inod 5). 



(3.3) 
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Table II. Pair-Correlation Function 



Regions 


Signs of (£0,4,4,4,^4) 








I & VI 


(+,+,+,-,-) & (-,-,-,+,+) 


9 


[14 


,14 


, 4 


, 4 


,|4|] 


II & VII 


(+,+,-,-,-) & (-,-,+,+,+) 




,14 


,14 


,14 


, I4|] 


III & VIII 


(+,+,-,-,+) k (-,-,+,+,-) 


^7[|4 


,14 


,141,14 


,|4|] 


IV & IX 


(+,-,-,-,+) & (-,+,+,+,-) 




,14 


,14 


,14 


, 141] 


V & X 


(+,-,-,+,+) & (-,+,+,-,-) 


9 


[14 


,14 


,14 


,14 


, 141] 



From Table [TTl we indeed find that this property holds. 



4. Enumeration of sites 

As mentioned in the introduction in the context of the Fibonacci Ising 
model, ''^'^ to calculate the q-dependent susceptibility (11 .ip for lattices for 
which the correlations are not translationally invariant, one needs to find 
a way to calculate suitable averages of the pair-correlation function. Since 
the meshes in the pentagrid — and even the distances between them — can be 
infinitesimally small in the thermodynamic limit, the problem of counting all 
the spin sites must be first solved. 

We proceed by considering in detail the parallelograms bounded by two 
sets of parallel grid lines in the pentagrid and examining all possible spin sites 
in each of these parallelograms. Let P{kj, kj^i) denote the parallelogram 
sandwiched between four grid lines kj — 1, kj, kj^i — 1 and kj+i for any j. 
(Throughout the entire paper, we let fcj+s = kj, i.e. the index j is considered 
mod 5.) Obviously, for all points z G P{kj, kj+i), we have Kj{z) = kj and 
Kj+i{z) = kj+i. The different choices of j give the different orientations of 
the parallelograms. Next, we determine how many spin sites a parallelogram 
may have, what are the integer vectors for these spins, etc. Since this section 
is rather lengthy, we have subdivided it into many parts, and put the main 
conclusion at the end. 
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4.1. Reference vector for P 

The vertices of parallelogram P{kj, fcj+i) can be calculated from (I2.3p as the 
intersections of grid lines in grids Gj and Gj+i, i.e. 



Gi n G.+i = <^ ^ G C 



^_ i[C^'(fc.+i-7,+i)-C^'+Hfc,-7,)] l 

sin(27r/5) J ' ^ ' ' 

for kj, kj^i G Z. Moreover, any point z in the interior of P{kj, kj^i) may be 
expressed in terms of e = (e^, ej+i), with ^ e^, ej+i ^ 1, as 

. = ^ (4.2) 

sin(27r/5) 

allowing us a change of notation Kj^^{e) = Kj^^{z{€)) for z E P. 

The four corners of parallelogram P{kj, kj^i) are given by e= (0, 0), (0, 1), 
(1, 0), or (1, 1) as can be seen from (14. ip . Now for each P{kj, kj+i), we pick 
a reference integer vector (fco, ■ ■ ■ , k^), which is related to the integer vector 
of the corner of P{kj, fcj+i) with e= (0, 0). Apart from the obvious identities 
kj = Kj{0) and fc^+i = Kj^i{0), we have 

k,+2 = i^,+2(0) 

= \P~\kj+i - 7j+i) - kj + 7i + 7j+2l = r«l - kj, 
kj+4 = Kj^i (0) 



\P \kj - -ij) - kj+i + 7^+1 + 7j-+4l = f/S] - kj+i, (4.3) 



in which 



a = a{kj+i) = p \kj+i - 7^+1) + 7^- + 7j-+2, 

/3 = Pikj) = p~\kj - 7,) + 7,-+i + 7,-+4. (4.4) 

However, for the last component of the reference integer vector we choose 

fc,+3 = 2 - H - [/?] =-[a\- [/3J ^ ir,+3(0), (4.5) 

where \_x\ denotes the "floor of x", which is the largest integer ^ x and 
[xj = \x~\ if and only if a; G Z. Since the pentagrid is regular, we find 
^ Z and the second equality in the above equation holds. The index of 
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any mesh, whose integer vector is K{z), is defined as Kj{z). It is shown 
by de Bruijn® that it has one of the four possible values, 1, 2, 3, or 4. We 
associate odd spins to meshes with index 1 or 3, and even spins to meshes 
with index 2 or 4. The index of the reference integer vector is 'Y^- kj = 2. 
From ([22D and (g^D we find 

Kj+;iO) = r-p-'(%+i - 7i+i + kj - 7j) + 7i+3)l = - P] ■ (4.6) 
Using 

{x} = x- [x\ , \-x-y] = - [x\ - [y\ + \-{x} - {y}] (4.7) 

and comparing (14.51) with (14. 6p . we find 

f A;,+3 - 1 for {a} + {p} ^ 1, 
[ for {a} + {(3}<1. 

This shows the mesh below the upper right corner with e = (0,0) belongs to 
the even sublattice and its integer vector is the reference vector of P only for 
{a} + {f3} < 1, but not for {a} + {f3} ^ 1. 

It would be more natural to choose this corner as our reference and to 
compare the integer vectors of other spins inside P = P{kj, kj+i) with it. 
This was what we did originally. However, we find that the rather odd 
choice of the reference vector given by (14.31) and (14. 5p . which may not even 
be the integer vector of a mesh, has made calculations much simpler. We 
next examine the differences between the integer vectors of the other spins 
in P{kj, kj^i) with respect to this reference vector. 



4.2. Integer vectors for z G P 

Substituting (14.21) into (12.41) and using (14.31) and (14.51) . we evaluate the integer 
vectors K{e) for every point in P that is not on a grid line of the pentagrid. 
We find that its components for m = 2, 3, 4 are given by 

Kj-f-^{e) = kjj^^ + (9-ft'j+in(e), 

dK,+^,{e) = \\+m{e) - 11 = LA,+m(e)J, (4.9) 

with 

Aj+2(e) = {a} + - p"^ej+i, Xj+ii.^) = {(3} + e^+i - p'^ej. (4.10) 
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and 

A,+3(e) = p-\e, + 6,+i) - {a} - {(3} + 1. (4.11) 

The last equality in (14.91) does not hold if Aj+m(e) is an integer, i.e. if the 
point is on a grid line of grid Gj+m. If we define the difference vector dK[e) 
for each mesh in P 

dK{e) = [dK^+2{e), dK,+s{e), dK,+,{e)] , (4.12) 

then for {a} + {(3} > 1, we have dK{0) = [0, -1, 0], and for {a} + < 1, 
dK{0) = [0,0,0]. 

For fixed 7j's with j = 0, ■ ■ ■ ,4, which are the shifts of the pentagrid, 
the d(/cj+i) and P{kj) in (14.40 are uniquely determined for each P{kj, kj+i). 
Consequently, the number of meshes in P and the difference vectors dK for 
each mesh are uniquely determined by (14.91) to (14.111) . The configurations of 
two parallelograms P and P' are the same, if they have the same number of 
meshes (spin sites), and the same sets of difference vectors. The difference in 
the configurations does not depend on the exact locations of the relevant grid 
lines or their intersections. However, whenever a grid line or an intersection 
moves in or out of the parallelogram P{kj, kj+i), the configuration changes. 

4.3. Relevant grid lines 

It is easy to see from (12.41) that dKj^^{e) changes its value whenever lines 
in the (j+m)th grid are crossed. Because ^ {x}, ej, e^+i < 1, we find from 
fICTD and fltnl) that 

- 1 < A,+2(e), A,+4(e) < 2, -1< A,+3(e) < 3. (4.13) 

Consequently, the only relevant grid lines for the parallelogram P{kj^kj+i) 
are those having integer labels — 1 + n', kj^^ — 1 + m, or ^^+4 — 1 + n, 
with n, n' = 0, 1 and m = 0, 1,2. Indeed, these are the only integer values 
that the Kj+m(e) in (14.91) . or /cj+m — 1 + r'^i+m(e)'| , can assume. The loci of 
these lines are given by linear equations in ej and e^+i as 

Aj+2(e) = n', Aj+3(e) = m, Xj+i{e) = n. (4.14) 

In our notation we suppress the two trivial components dKj{e) = ^Kj-^-l(e) = 0, since 
in parallelogram P we have by definition Kj{e) = kj, Kj^i{e) = fcj+i. 
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From (14.101) . we find 

< Aj+2(e) < 2 if {a} > p-'\ 

< Xj+4{e) < 2 if {(3} > p-\ (4.15) 

Therefore the equation Aj+2 = (Aj+4 = 0) cannot be satisfied if {a} > 
{{P} > p^^), while equations = 1 and Aj+4 = 1 always have solutions 
in P. This means that the only grid lines in grids Gj+2 and Gj+4 crossing P 
are given by 

kj^2, G P always, 
/cj+s-IgP if {a}<p'^, 
kj+i-leP if {P}<p-\ (4.16) 

From (14. lip we find that the grid lines /c^+s + m are parallel to the diagonal 
ej + ej+i = 1, and 

l-p-3 < A,+3(e) < 2p-i + l if < {a} + {13} < p-^ 

< A,+3(e) < 2 if p-^ < {a} + {(3} < 1, 

< A,+3(e) < 2p-' if 1< {«} + m < 2p~\ 

-1 < A,+3(e) < 1 if 2p-^ < {a} + {f3} < 2, (4.17) 

where we used the identity p~^ = 2p~^ — 1. This means, there can be at most 
two lines of grid Gj^^ going through the inside of the parallelogram. We find 
that Aj_|_3 = 2 has a solution in P if < {a} + {/?} < Aj+s = 1 has a 
solution in P for < {a} + {[3} < 2p~^; Aj+3 = has a solution in P for 
1 < {a} + {j3}. These facts can be summarized as 

v^-2. 



kj+3 + 1 e P 


if 


{a} + m < p-' = 


kj+3 e P 


if 


{a} + m < 2p-' 


kj+3 -leP 


if 


{a} + {(3}>l. 



(4.18) 

At this point one may note the symmetry under kj^2 ^ kj^4^ and a <-> /5 in 
conditions and (KIM . 
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4.4. Intersections of the grid lines 

Next, we need to calculate the positions of the intersections. We let a"'™ 
denote the intersection of a pair of grid lines in G^+a and Gj+4 numbered 
kj^-^ — l+m and kj^4 — l+n, 6"'"^ the intersection of lines kj^^ — l + m and 
fcj+2 — 1+^' in Gj+3 and Gj+2, while c"'" the intersection of lines 1+n 
and — ^ + in Gj+4 and Gj+2- The locations of these intersections are 
found by solving the corresponding linear equations given in f l4.14p . We find 

a^,ni^ (p({/3}_^)+m+^_l + {a},p-i({a} + n+m-l)), (4.19) 
_ + p^{a}-n')+m+n'-l + {(3}), (4.20) 

c"'"' = (pn'+n-j9{a}-{/5}, pn+n' -{a}-p{(3}). (4.21) 

Clearly, whenever both components (e-,-, e^+i) of an intersection are positive 
and less than 1, it is inside the parallelogram P{kj, kj+i). This way, we find 
from fl4.19p the conditions for the three possible cases, namely 

aO'i gP(A;„A;,+i) if {[3} < p-' ~ p-' {a} , 

a^'° G P{kj, kj+i) if {/?} > 1 - p'^{a}, 

f < |a| < p^^ and 
a^'i G P(L-,A;.+i) if <^ (4.22) 
' ' \ p-'-p-'{a}<{P}<l-p-'{a}. 

Similarly, from (14.201) we obtain 

6°'^ gP(%,A;,+i) if {[3}<l-p{a}, 

b'-"^ eP{k„k,+,) if {P}>p{l-{a}), 

{ < {/?} < p-^ and 
b^'^GP(A:,,A:,+0 if ^ (4.23) 

and from (I4.2ip we get 

cO'i gP(A;„A;,+i) if p-' -p{a} <{f3} <p-\l-{a}), 

c^'° e P{kj,k^+i) if p-2-p-H"}<{/5}<l-J5{a}, 

c^'i G P(A;j, fcj+i) if {p} > max{l~p-^{a},p~p{a}). (4.24) 
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Note the symmetry between (14.221) and fl4.23p and between the first two hnes 
of fl4.24p under a ^ (3 and implicitly kj+2 ^ 

Three lines cannot exactly meet in a common intersection, as was shown 
by de Bruijn® for a regular pentagrid. However, they can meet arbitrarily 
close and the theoretical limiting conditions of triple intersection 







= ci'i " 















play a role in the following subsection. 



{«} + {/?} =p, 

{a} + {P}=p-\ 

{a} + {P}=p-\ (4.25) 



4.5. The twenty-four allowed configurations 

Next, we use ( KW) . fl4A8ll . (141221) - ( 14:251) to study how the configuration 
C{m) of parallelogram P{kj, fcj+i) depends on the values of {a} = {d(/cj+i)} 
and {/?} = {P{kj)}. We show the various cases in Fig. E]for j = 0. For jj^O 
we need to rotate each picture j times 72°. 

For {«},{/?} > p~^, the three grid lines %+2, %+3 — 1, and their 
three intersections a^'°, b^'*^, c^'^ are inside the parallelogram, producing 
seven meshes in P{kj, A^j+i), as shown in Figs. [6](a) and[6](b). The difference 
between the two cases is that the intersection c^'^ is on opposite sides of the 
grid line 1, as the sign of {a} + {P}—p changes, cf. (I4.25p . The grid line 
kj+3 — l moves upward toward the upper right corner as {/?} decreases. It is 
below the diagonal ej+Cj+i = 1 for + > p, corresponding to Fig.E](a), 
and above the diagonal for {a} + {/?} < p, as shown in Fig. [6](b). The index 
of the inner triangle changes from odd to even in view of (I4.10p . (14. lip and 
(14.140 . Hence, the spin configurations for the two cases are different: C(l) in 
Fig. [n](a) has 4 odd sites and 3 even sites, and C(2) shown in Fig. E](b) has 
3 odd sites and 4 even sites. Also, C(l) is the only configuration for which 
the reference integer vector does not correspond to an actual mesh. 

For {a} > p~^ and < {j3} < p~^, we have six cases denoted by C(3) 
to C(8) arranged in the decreasing order of {/?}. For 1 — p~^{a} < {/?}, 
the configuration C(3) is almost the same as C(2) except having one more 
even site because line kj+4 — l is now inside the parallelogram, as shown in 
Fig. [6](c). As {P} decreases, line fcj_|_4 moves downward and for 2p~^ — {a} < 
{P} < 1 — p~^{a}, the intersections of line with lines fcj+a — 1 (a^'°) and 
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Fig. 6. The spin configurations C(m), m — 1,...,24, for parallelogram P{ko,ki) are 
shown in (a) through (x), where N denotes the total number of spin sites inside P{ko, ki), 
No the number of odd spins sites, and the number of even spin sites. The mesh that 
contains the dot is the special even site whose integer vector is the reference integer vector 
for P(fco, fci), such that dK = (0, 0, 0). For {a} + {/?} > p, such a site does not exist, as 
seen in 6(a). For {a} + {/?} < 1, this site is the mesh right below the upper right corner 
of P. 
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Fig. 6. Continued: The remaining twelve configurations. The number of spin sites N 
varies between 6 and 12, whereas and iVo vary between 3 and 7. 
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kj+2 (c^'^) are seen from (14.221) and (14 .24^ to be outside of P{kj, fcj+i). As a 
result, two of the even sites are now outside, and C(4), shown in Fig. [6](d), 
has only six sites, three of which are odd, and three even. 

As {P} decreases further top(l— {a}) < {P} < 2p~^—{a}, line /c^+s is now 
inside P{kj, kj^i), as seen from (I4.18p . Thus, configuration C(5), shown in 
Fig. E](e), has one more even site than C(4). Both lines and /cj+3— 1 move 
upward as {P} decreases. For 1 — {a} < < p{l — {a}), the intersection 
b^'^ of line kj+3 with kj-^.2 moves inside of P{kj, /cj+i), while 6^'°, which is the 
intersection of 1 and moves out, giving rise to configuration C(6) 

shown in Fig. [6](f), with three even sites and four odd sites. 

For {l—{a}) < {/?} < 1— {a}, line kj+3—l moves out of P, as shown in 
Fig.[6](g). Its configuration C(7) has 3 odd sites and 3 even sites differing from 
C(6) in that the odd site with dK{0) = (0, —1, 0) is now outside P{kj, kj+i). 
For < {P} < p~^{l—{a}), the intersections a*^'^ and c°'^ of line kj+4 — l with 
kj+s and are now seen from (I4.22p and (14.241) to be inside of P{kj, kj-^i) 
adding 2 more odd sites to C(8), which is shown in Fig. [H](h). 

In Fig. [6](i) through Fig. [6](n), the six cases C(9) through C(14) are 
shown for {a} < p^^ and {P} > p~^. Because of (14.31) and (14.41) we can use 
the reflection symmetry {a} <-> kj ^ fc^+i, ^ kj^^, which was 
noted also in the previous subsections. Thus these cases are similar to the 
configurations C(3) to C(8), and obtainable simply by replacing kj^^ + n by 
kj+2 + n and vice versa. To summarize, we find C(9) has 3 odd sites and 5 
even sites; C(10) 3 odd sites and 3 even sites; C(ll) 3 odd sites and 4 even 
sites; C(12) 4 odd sites and 3 even sites; C(13) has 3 odd sites and 3 even 
sites, and C(14) 5 odd sites and 3 even sites. 

For {a}, {P} < p~^, at least five grid lines fcj+2 — 1, kj+4 and 

kj+4—1 are inside P{kj, kj+i). In Fig.[6](o), we show configuration C(15) valid 
for 1 < {q;} + {/5}, when both lines — 1 and fcj+3 and the intersections b^'^ 
and a^'^ are inside P{kj, fcj+i), as seen from (I4.22p and (I4.23p . Configuration 
C(15) has 4 odd sites and 5 even sites. 

For the region satisfying the three inequalities {P} < ! — {«}, {/?} > 
l—p{a} and {P} > p~^{l — {a}), grid line /cj+3 and the corresponding odd 
site with dK = (0, —1, 0) are now outside P{kj, fcj+i), such that configuration 
C(16) shown in Fig. [H](p) has one site less than C(15). It has 3 odd and 5 
even sites. 

For l—p{a} < {P} < p~^{l — {a}), the intersections a°'^ and c°'^ of 
line — 1 are both also inside P{kj, fcj+i), as seen from (I4.22p and (14.241) . 
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adding two odd sites to C(16). Thus, C(17) in Fig. [6](q) has 5 odd sites and 
5 even. However, for p^^{l—{a}) < {/3} < l—p{a} the intersections 6°'^ and 
c^'° of hnes — 1 are now inside instead, adding two different odd sites to 
configuration C(16). The resulting configuration C(18) shown in Fig. [6](r) 
also has ten sites and relates to C(17) by the above reflection symmetry. 

When both conditions p^^—p{a} < < l—p{a} and p^"^ —p^^{a} < 
{P} < p^^{l—{a}) are satisfied, we find six intersections a^'^, b^'^, a^'^, b^''^, 
c^'^ and c^'" inside P{kj, kj+i). As a result, there are 12 sites for the two 
cases C(19) and C(20) shown in Fig. [6](s) and Fig. [6](t) respectively. The 
difference between the two cases is that the intersections c^'^ and c^'° are on 
opposite sides of the grid line as the sign of {a} + {/?} — p~^ changes, 
cf. f l4.25p . For {a} + {P} > p~^, line fcj+a lies below the diagonal, as can be 
seen from (14.111) and (14.141) . so that its configuration C(19) has 7 odd sites 
and 5 even sites; for {a} + {P} < p~^ , fcj+s is above the diagonal, and C(20) 
has 5 odd sites and 7 even sites. 

For p~^ —p{a} < {P} < p~'^ — p~^ {<y} , intersections a^'^ and c^'° are 
no longer inside P{kj, kj+i). As a consequence two of the even sites are 
now outside, leaving C(21) shown in Fig. [H](u) with 5 odd sites and 5 even. 
For p~'^—p~^{a} < {P} < p^^—p{oi}, however, intersections 6^'^ and c°'^ 
are outside P{kj, kj+i) instead, such that two different even sites are now 
outside. The resulting configuration C(22) shown in Fig. [6](v) has the same 
5 odd sites as C(21), but a different set of 5 even sites. Again, C(21) and 
C(22) are related by the aforementioned reflection symmetry. 

For {P} < p^'^—p^^{a} and {/?} < p~^—p{a}, only two intersections 
a°'^ and b^'^ are still inside P{kj, fcj+i). These are the cases C(23) and C(24) 
shown in Fig. [6](w) and Fig. [6](x). For {P} + {a} > p~^, configuration 
C(23), shown in Fig. [ni(w), has eight sites: 5 odd and 3 even. Finally, for 
< {P} < p~^ — {a}, the grid line + 1 is inside P, so that C(24) shown 
in Fig. [6](x) has nine sites, of which the 5 odd sites are identical to those of 
C(20) through C(23). 

In summary, the boundaries for the above 24 regions are the 13 lines given 



Q-Dependent Susceptibility on Penrose Tiles 



25 



by 

{a}=p-\ m=v-\ 

+ {a} = p-i, 1, 2p-\ orp, 
{(3}+p-^{a}=p-\ p-\ or 1, 

{(3} + p{a} = p-^ , 1, orp, (4.26) 

which are also the boundaries of the inequahties in f l4.16p . (14.181) . (14.221) . 
(14.231) and (14.241) . together with the two conditions in (14.251) . These are 
exactly all conditions for three grid lines to meet at a corner of P, on an edge 
of P, or inside of P, respectively, i.e. the only conditions under which some 
mesh can appear or disappear in P under shifts of grid lines. 

In Fig. [71(a), we plot the boundaries lines given by (14.261) in the unit square 
with {a} and {/?} along the horizontal and vertical axes. These lines indeed 
divide the unit square into 24 regions. Each of these 24 regions corresponds to 
a different configuration of the parallelogram P(/cj, kj+i). The above analysis 
shows that the parallelograms can only have 6, 7, 8, 9, 10 or 12 sites inside. 
The position of {d(A;j+i)} and {l3{kj)} in the unit square shown in Fig. [7](a) 
determines which configuration the parallelog ram P(^kj, kj^ij is in. 

The areas of the 24 regions in Fig. [7] (a) can be easily calculated. We find, 
after using the formula of the area of a triangle in terms of the coordinates 
of its three vertices, 

A{1) = A{2) = I-Ip= \p-\ 

A(3) = A{b) = A{6) = A{8) = A{9) = A{11) = A{12) 

= A{U) = A{19) = A{20) = |p - 4 = 
A{4) = A{7) = A{10) = A{13) = A{15) = A{17) = A{18) 

= A{21) = A{22) = A{2A) = f - Ap = \p~\ 
A{IQ) = A{23) = 9p-f = - (4.27) 

The areas of all 22 triangular areas differ by powers of the golden ratio p. 

For later reference, we display in Fig. [7](b) the eight different regions 
with equivalent odd configurations. Two odd configurations are equivalent. 
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Fig. 7. (a) The 24 regions for the 24 different configurations are shown within the unit 
square with {a} and {/?} along the horizontal and vertical axes. The integer m denotes 
the mth region for the mth configuration, (b) The 8 regions for the 8 different odd 
configurations are shown. 



if they have an equal number of odd sites with the same sets of integer vectors 
for these odd sites. This identification can most easily be made using the 
information in Tables 111(a) and 111(b) below. 



4.6. Probability 

From the definition (14.41) . we find that the fractional parts {a{kj^i)} and 
{P{kj)} are related to the golden ratio p which is irrational. From the well- 
known theorem of Kronecker,'E3' we conclude that {d(/cj+i)} and {P{kj)} 
are everywhere dense and uniformly distributed in the interval (0, 1), as the 
integers kj and kj+i vary from — oo to oo. As a consequence, every point 
in the unit square in Fig. [7] is equally probable. Therefore, the frequency or 
probability for a parallelogram to be in one of the twenty-four configurations, 
say m, is given by the area A(m) of the mth region. 

Although the pentagrids are different for different choices of the shifts jj, 
and the values of {^(fcj+i)} and {P{kj)} are also different for different 7j's, 
this does not change the probability distributions of {a{kj+i)} and {P{kj)} 
in the thermodynamic limit of kj and kj^i varying from — oo to oo. In other 
words, the area A{m) for the mth configuration is independent of the j/s, 
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and is the same for all regular pentagrids. 



4.7. Difference vectors 



For {a} + {/?} > 1 the grid line kjj^^ — l is seen from fl4.18p to be inside 
P{kj,kj+i). Consequently, the mesh below corner e = (0,0) is odd and 
dK = (0,-1,0). For {a} + < 1, line /cj+3 — 1 is outside P{kj,kj+i), 
and the mesh at corner e = (0, 0) is even. In this case, its integer vector is 
identical to the reference vector for P{kj,kj+i) with dK{0) = (0,0,0). The 
difference vectors dK for all other meshes in the parallelograms in each of 
the twenty-four configurations can be easily obtained from Fig. El as dKj+ja 
changes its values only when a lines in the (j + m)th grid is crossed. 



In Tables Ill(a) and Ill(b) , we list for each of the 24 configurations all 
the difference vectors dK = {dKj+2,dKj+3,dKj^4), with Table Ill(a) for 



{a} + {(3} > 1 and Table Ill(b) for {a} + {(3} < 1. From these tables one 
can also immediately read off which regions have equivalent odd (or even) 
configurations. 



4.8. Average Number A/" 

We shall now calculate A/", which is the average number of spin sites in a 
parallelogram. The total number of sites can be evaluated by counting all 
the sites in each parallelograms P{kj, fcj+i), and then adding all of them to- 
gether for all the P's. This is equivalent to splitting the summation over all 
sites into two parts — first summing over all sites in P represented by their 
integer vectors K{e) and then adding all of them for all the parallelograms. 
Let there be Ai lines in each of the five grids, so that there are parallel- 
ograms, ignoring boundary effects that cancel in the thermodynamic limit. 
The average J\f then equals the total number of lattice sites divided by A^^. 

We have already shown that each parallelogram P{kj, kj^i) is in one of 
24 configurations C(m) uniquely determined by the values of {a{kj+i)} and 
{P{kj)}. The allowed configurations C(m) have N{ra) = 6, 7, 8, 9, 10 or 12 
sites inside P. The frequency or probability A{m) is defined as the number 
of parallelograms in the mth configuration divided by the total number of 
parallelograms. If we let kj and fcj+i in P{kj,kj+i) run over the A4 values, 
each of the parallelograms is counted once. 

As Ai approaches 00, so that —00 < kj, kj+i < 00, the values of {^(fcj+i)} 
and {${kj)} are everywhere dense and uniformly distributed® between 



Table 1 1 1(a). {ct} + {(3} > 1 





C(l) 


C(2) 


C(3) 


C(4) 


C(5) 


C(6) 


C(9) 


C(10) 


C(ll) 


C(12) 


C(15) 




[0,-1,0] 


[0,-1,0] 


[0,-1,0] 


[0,-1,0] 


[0,-1,0] 


[0,-1,0] 


[0,-1,0] 


[0,-1,0] 


[0,-1,0] 


[0,-1,0] 


[0,-1,0] 


o 


[1,0,0] 


[1,0,0] 


[1,0,0] 


[1,0,0] 


[1,0,0] 


[1,0,0] 


[1,0,0] 


[1,0,0] 


[1,0, 0[ 


[1,0,0] 


[1,0,0] 


d 


[0,0.1] 


[0.0,1] 


[0.0,1] 


[0.0,1] 


[0.0,1] 


[0,0,1] 


]0.0,1] 


[0,0.1] 


[0,0, 1[ 


[0.0,1] 


[0.0,1] 


d 


[1-1-1] 










[0,1,0] 








[0.1,0] 


[0.1,0] 




[1,-1,0] 


[1,-1,0] 


[1,-1,0] 


[1,-1,0] 


[1,-1,0] 




[1,-1,0] 










e 




[0, 0, 0] 


[0, 0, 0] 


[0, 0, 0] 


[0,0,0] 


[0,0,0] 


[0,0,0] 


[0,0,0] 


[0,0, 0[ 


[0, 0, 0] 


[0, 0, 0] 


V 


[O.-l.ll 


[0,-1,1] 


[0,-1,1] 








[0,-1,1] 


[0,-1,1] 


[O.-l. 1] 






e 


[1,0.11 


[1.0,1] 


[1.0,1] 








[1.0,1] 










n 






[1,0.-1] 


[1,0.-1] 


[1,0.-1] 
]1.1,0] 


[1,0,-1] 
[1,1,0] 


[-1,0,1] 


[-1,0,11 


[-1,0,1] 
[0,1,1] 


[-1,0,1] 
[0,1,1] 


[1,0,-1] 
]-l,0,l] 
[1,1,0] 
[0,1,1] 



Table lll(b). {a} + {(3} < 1 





C(7) 


C(8) 


c(i;i) 


C(14) 


C(16) 


C(17) 


C(18) 


C(19) 


C(20) 


C(21) 


C(22) 


C(2.'i) 


C(24) 




[1,0, 0| 


[1,0,0] 


[1,0,0] 


[1,0, 0| 


[1,0, 0| 


[1,0,0] 


[1,0, 0| 


[1,0, 0| 














[0,1, 0| 


[0.1,0] 


[0,1,0] 


[0,1, 0| 


[0,1.01 


[0,1,0] 


[0,1, 0| 


[0,1. 0| 


[0,1,0] 


[0.1,0] 


[0,1.01 


[0,1, 0| 


[0.1,0] 


o 


[0, 0. 1| 


[0. 0, 1] 


[0. 0, 1] 


[0, 0. 1| 


[0,0. 1| 


[0. 0, 1] 


[0, 0, 1| 


[0, 0. 1| 












d 




[0,0.-1] 








[0.0, -1| 




[0,0.-1] 


[0,0.-1] 


[0.0, -1| 


[0.0,-1] 


[0,0.-1] 


[0.0, -1| 


d 




[1,1,-1] 








[l.l,-ll 




[1,1.-1] 


[1,1.-1] 


[l.l.-l] 


[1.1,-1] 


[1,1.-1] 


[l.l,-ll 










[-1,0,01 






[-1,0,01 


[-1.0,01 


[-1,0,0] 


[-1,0,01 


[-1,0,01 


[-1,0,0] 


[-1,0,0] 










[-1.1.1] 






1-1,1,1] 


[-1.1,1] 


[-1.1,1] 


1-1,1,1] 


1-1,1,1] 


[-1,1,1] 


1-1,1,1] 




[0,0,0] 


[0,0,0] 


[0,0,0] 


[0,0,0] 


[0,0,0] 


[0,0,0] 


[0,0,0] 


[0, 0, 0] 


[0,0,0] 


[0,0,0] 


[0,0, 0| 


[0,0,0] 


10, 0, 0] 








[0.1,1] 


[0,1. 1| 


[0,1. 1| 


[0.1,1] 


[0,1, 1| 


[0,1. 1| 


[0,1,1] 




[0,1, 1| 






e 


[1,1,0] 


[1.1,0] 






[1,1.01 


[1,1.0] 


[1,1,01 


[1,1.01 


[1.1,0] 


[1.1,0] 








V 


















[-1,1,0] 


[-1,1,01 


[-1,1,0] 


[-1.1,0] 


[-1,1,0] 


e 


















[0,1,-1] 


10,1,-1] 


[0,1,-1] 


[0,1,-1] 


[0,1,-1] 


n 


[1,0,-1] 


[1,0,-1] 






[1,0,-1] 


[1,0,-1] 


[1,0,-1] 


[1,0,-1] 


[1,0,-1] 


11,0,-1] 














[-1,0,11 


[-1,0,11 


[-1,0,1] 


[-1,0,1] 


[-1,0,11 


[-1,0,11 


[-1,0,1] 




1-1,0,1] 
































10, 2, 0] 
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and 1. The frequency A{m) is the area of the mth region in the unit square 
in Fig. [71(a). The values of these A(m)'s are hsted in (14.271) . Denoting the 
number of sites in the mth configuration by N{m), with values also given in 
the captions of Fig. [6l then the average number of sites per parallelogram is 

K{z£C) allP K{zGP) 



M- 



24 N{m) 24 

A{m) XI ^ = XI Mm)N{m) = 5p. (4.29) 

m=l n=l m=l 



where we have used the notation K{z G C) to denote the integer vectors of 
all the meshes in the pentagrid, while K{z G P) denotes only those meshes 
in parallelogram P. 

4.9. Penrose Tiles 

The above method provides an alternative way to draw the Penrose tiles. To 
illustrate this, we let j = 0. For some fixed shifts 7j, we let — J ^ /eq, ki ^ J 
for some positive integer J. For any values /cq and ki in this set, {a{ki)}, 
{P{ko)} are uniquely determined from (14.41) . The elements k^, for m = 2, 
3 ,4, of the reference vector in P{ko, ki) are given by (14. 3 p and (14.51) . From 
the values of {a{ki)} and {/3(fco)}, we can determine from Fig. [7](a), what 
configuration C{m) the parallelogram P{ko, ki) is in. Then we can use Tables 



Ill(a) and Ill(b) to obtain the difference vectors dK for the N{m) sites inside 
P{ko,ki). The integer vectors for these different sites in P{ko,ki) are then 
given by 

i?(e) = {ko, fci, k2 + dK2, ks + dK^, k^ + dK^). (4.30) 

We use (12.51) and (I4.30p to obtain the positions of the spins in the complex 
plane for the meshes in P(fco, fci). Hence, as /cq and ki run over the values 
from — J to J we obtain the positions of the spins in both odd and even 
sublattices shown in Fig. [5l This figure has been plotted using Maple. 

4.10. Summary 

Consider the parallelograms P{kj, kj+i) which contains all points ^ G C 
such that Kj{z) = kj and fC,+i(2;) = fcj+i, cf. (12. 4p . The configurations 
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of two such parallelograms are considered to be the same, if they contain 
the same number of spin sites and the corresponding sites have the same 
difference vectors. The different configurations do not depend on the exact 
locations of the relevant grid lines or their intersections. However, whenever 
a grid line or an intersection moves in or out of the parallelogram P, the 
configuration changes. The above analysis shows that there are only 24 
allowed configurations, with 6, 7, 8, 9, 10 and 12 sites inside P. 

The configuration of kj+i) is uniquely determined by the values of 
{a{kj^i)} and defined in f l4.4p . By examining the locations of the 

relevant grid lines and their intersections, we find that the unit square, with 
{a} and {(3} along the horizontal and vertical axes, is divided into 24 regions, 
corresponding to the 24 possible configurations of the parallelogram P. This 
is shown Fig. [7] (a), and the 24 configurations of P are shown in Fig. O 
Using the theorem of Kronecker,''^ we find that the area A{m) of the mth 
region is actually proportional to the probability for the mth configuration to 
occur. Even though the pentagrids and configurations of the A;j+i)'s are 
different for different choices of the shifts 7j, the area A{m) is independent 
of these shifts and is the same for all regular pentagrids. 

In each parallelogram P{kj, kj^i), a reference integer vector is chosen 
whose components are given by (14. 3 p and (14. 5p . The difference vectors with 
respect to this reference vector are defined in (14.90 and calculated for all 
of the sites inside P. The number of sites N{m), the area A{m), and the 
difference vectors dK for the N{m) spin sites for m = 1...24 are listed in 
Fig. El Eq. Km and Tables [in(t)] and |III(b)| respectively. 



5. Susceptibility 

There are three Z-invariant Ising models that can be defined on the vertices 
of the Penrose rhombus tiles using the prescriptions of Section [2l Model 1 
has spins on all odd sites only, interacting along the diagonals of the tiles, 
as is illustrated in Fig. [51 Model 2 is defined similarly with spins only on 
the even sites. Model 3 has all sites of the Penrose tiling, but the even and 
odd sites are decoupled, with the odd spins interacting as in model 1 and the 
even spins as in model 2. We will see that the three models have the identical 
wavevector-dependent susceptibility P^r spin site in the thermodynamic 
limit. 

The physical positions of the spins have been expressed in (12.51) as complex 
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numbers depending on the integer vectors K{z) of the meshes. Let g be a 
complex number and q = Qx + iQy so that q* denotes its complex conjugate 
(while q = {qx,qy)), then the q-dependent susceptibility is 

A:BTx(q)=J,im^-^ 5^ E cosRejg* - } 

K(zeC) K{z'eC) i=o 

2x°(q), (model 1), 

X [{('KizfKiz')) - {^Kiz)){^Kiz'))] = { 2x'(q), (model 2), (5.1) 

X°(q) + X''(q), (models). 

Here the double sums denoted by K{z G C) are over all the odd spin sites in 
model 1, over all the even sites in model 2, and over all sites in model 3. In 
the last case, since the spins on the odd sublattice do not interact with those 
on the even sublattice, the q-dependent susceptibility x(q) becomes the sum 
of two parts: X°(q) denoting the contribution from the odd sublattice and 
X''(q) from the even sublattice. For model 3 the average number of spin sites 
per parallelogram was given in fl4.29p as A/" = 5p. For models 1 and 2 this 
number becomes A/" = 5p/2, explaining the extra factors 2 in the last member 
of (lO). 

We shall first consider x°{q) and show later that x'^(q) is equal to it, 
implying that the susceptibilities of the three models are indeed equal. 

5.1. Calculation of x°(q) 

We again split the sum over all odd spin sites into two parts as in (14.281) . and 
let K°{z G P) run over all odd spins in parallelogram P. Consequently, 

*°<i)=i°Lm<5 5: E E E m-(z%K%z)), (5,2) 

allP K°{zGP) all P' K°{z'£P') 

where 

4 

X cosRe[q*J2[Ki^')-Ki^)]C}^ (5-3) 



n=0 
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with {(yiy'Y denoting the connected pair correlation function, subtracting the 
contribution from the spontaneous magnetization. 

Now we let P = P{kj,kj+i) and P' = P{kj + i,kj+i + i'). As ^ cx), 
£ and i' are kept fixed, and kj and kj^i vary from — cxo to oo, all the P's 
and P"s are counted once. It is also evident that the different choices of j 
correspond to choosing one of five orientations for the parallelograms, and 
they should give the same q-dependent susceptibility. 

In (14.21) we defined z = z{e) for z & P. Similarly, for z' G P', we let 

so that z' ^ e' = (e^-, e^_j^^) and < ej,ej+i < 1- The corresponding integer 

vectors K{z') = K{z'{e')) = K'{e') also have two fixed components Kj{e') = 
kj+i = k'j and K'j^^{e') = = k'j^^. Now, following (lOIl . we let 

a' = a{k'j^^) = p^^{kj+i+i' --fj+i) + 7^- + 7^+2 = a + p~^i', 

(3' ^ Pik'^) = p-\k,+i-^,) + 7,+4 + 7,+i = P + p-'t (5.5) 

According to (14.31) and (14. 5 p the reference integer vector k' for P' is chosen 
to have components, 

-+4 = - k,+i -i' = A:,+4 + 5,+4 + [p-'i\ - i', (5.6) 

in whicljfl 

= [{"} + {P''n\ , ^.+4 = m + {p-'i}\ . (5.7) 
Substituting (EJD into ( 12^ we obtain for m = 2, 3, 4 

K'^U^') = k'^^^ + dK'^^^{e'), dK'^_,Je') = [X'^U^')\, (5.8) 



8 From the first members of ([5^ and (gSl) we find that Sj+2 = [a'] - \a\ - [p'^t'l = 
\a+p-^e']-\a\-lp-H'\ = {{aj-l + ip-^l'jl = [{a} + {p-^f }J . For the last two steps 
here and the remaining steps in the derivation of (|5.6p and (|5.7|) we must make explicit 
use of the fact that a, a' , P and f3' are not integer for a regular pentagrid. 
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where 



a;+2(^') = {«'} + 4-^ 4+1' a;+4(^') = {/?'} + 4+1 4- (5-9) 

and 

a;+3(^') = p-\^', + s+i) - {«'} - if^'} + 1- (5-10) 

Comparing (15.81) to (15.101) with (14. 9 p to (14. lip , we find that the dependence 
of A^^_jn(e') on {a'} and {/?'} is the same as the dependence of Aj+m(e) on {a} 
and Consequently, the configurations of P' depend on {a'} and {j3'} 

in the same way as the configurations of P on {a} and {/?}. Therefore, the 
position of {a'} and {/?'} in the unit square in Fig. [7](a) with {«'} and {j3'} 
along the horizontal and vertical axes uniquely determines the configuration 
of P'. The difference vectors dK'(e') for the sites in P', which is in one of 



the 24 configurations, are again given in Tables Ill(a) and Ill(b) 



The above results are valid for all spin configurations in P and P', but 



we shall consider only the odd spins at first. By examining Tables Ill(a) and 



Ill(b) , we can find that the odd spin configurations of the parallelograms 
are simpler, because several connected regions — see C(2) to C(5), or C(9) 
to C(ll) in Fig. [7] (a) as examples — have the same odd spin configurations. 
In fact, there are only eight distinct odd spin configurations. In Fig. [7](b), 
the regions for these 8 odd spin configurations are shown in the unit square 
whose axes are the {a} and {(3} directions. 

Listed in Table [IV] are the number of odd sites N{m) in P, the region 
of validity R{m), and the area A{m) of the mth odd configuration for all 
m = 1, ■ ■ ■ , 8. In the mth odd configuration, the difference vectors of the 
N{m) spins are denoted by dK^"^'^^ for n = 1, - ■ ■ , N{m). They are equal to 
the difference vectors dK'^ of the odd sites in some configuration C(/), listed 



in Tables Ill(a) and Ill(b) for / = 1, ■ ■ ■ , 24. In the last column of Table IIVI 
it is indicated which C(/)'s correspond to a given m. 

Let the distances £ and £' between the two parallelograms P and P' be 
fixed, but kj and kj+i vary from — oo to oo. Then {a'} and {P'} given by (15. 5p 
are also everj^vhere dense and uniformly distributed in the interval (0,1). 
The area A(m) is again the probability or frequency of the mth configuration. 



Table IV. The eight regions for the odd sublattice 



m 


N{m) 


A{m) 


R{m) 


1 ^ n ^ N{m) 


1 


4 




p-^ <{a}<lkp- {a} <{/3} <1 


dK'' in C(l) 


2 


3 




< {a} ^ p-i & 1 - p-^a] < {fn < 1; 

p-^ ^{a} <lk p{l - {a}) < {P} <p - {a} 


dK° in C(2) or 

V / 

C(/), /=3,4,5,9, 10,11 


3 


4 




< |a| ^ p-^ & 1 - \a} <\p} <1- p~^\a}: 
^ {«} < 1 & 1 - {a} < {f3} < p{l - {a}) 


in C(6), 
C(12) or C(15) 


4 


3 




< {«} ^ & 1 - p{a} < {/?} < 1 - {a}; 
^ {a} < 1 & p-\l - {a}) <{f3}<l- {a} 


dK° in C(7), 
C(13) or C(16) 


5 


5 


1 A 


0<{P}< p-^ & p-\l - {13}) <{«}<!- p{(5} 


in C(8) or C(17) 


6 


5 




< {a} < & - {a}) <{/?}<! -p{a} 


inC(14) orC(18) 


7 


7 




< {a} ^ p-2 & p-i - {a} < {/?} < p-\l - {a}); 
p"^ ^ {a} < p~^ & - {a} < {/?} < 1 - p{a} 


dK° in C(19) 


8 


5 




< {a} < p-^ & < {/?} < - {a} 


in C(20) to C(24) 
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From fl5.5p we find 

{a} + a for {a} + a < 1 
{a} + a — 1 for {a} + a ^ 1 

{/?} + b for + 6 < 1 
{l3} + b-l for {/?} + 6 ^ 1 

The plot of the eight regions for the odd configurations of P = P{kj, kj^i 



{a'} 

m 



b={p-H}. (5.11) 



a 



{PI 



^2 






5^ 






1 




\ 1 




6 \ \\ 
V \ \ 3 


2 


2 


8 






{a} 

(a) 

Fig. 8. (a) The eight odd spin configurations in P(kj+ 2, plotted with {a} 

and {/?} as the horizontal and vertical axes, (b) Overlapping the plot in part (a) with 
Fig. 7 (b) we can determine the joint probability Am^m'i'^, 3) for P{kj, fcj+i) to be in the 
mth configuration and P(fcj + 2, fcj+i +3) in m'th configuration, with m,m' — 1, • • • ,8, 
geometrically. The area of the shaded region represents the probability for m = 8 and 
m' = 2. 

has been given in Fig. [7](b). The plot for P' = P{kj+i, kj-^i+i') is the same, 
only now with {a'} and {(3'} along the axes. However, if the eight regions are 
replotted with {a} and {/?} along the horizontal and vertical axes, then we 
obtain unit squares as shown in Fig. [8] (a) [for the special case of £ = 2 and 
i' = 3]. Because of the relation (15. lip , we find that Fig. [8](a) can be obtained 
from Fig. [7](b) by cutting a horizontal slice with width b from the bottom 
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of Fig. [71(b), and pasting it on the top; then cutting a vertical shce of width 
a from the left and pasting it to the right. After the cutting and pasting, 
the connected region R{m) for P' in Fig. [7](b) becomes R'{m) in Fig. [8](a), 
which — for £ 7^ or 7^ — may consist of disjointed pieces pasted in up 
to four different sections of the unit square. The probability for P' to be in 
the mth configuration is still the area A(m) = area(i?(m)) = area(i?'(m)), 
which in the latter case, could be a sum of areas of disjointed pieces. 

For fixed i and I' [chosen to be £ = 2 and ^' = 3 in Fig. [8](a)], the 
position of a = and /? = l^ikj) in the unit square shown in Fig. [71(b) 

uniquely determines the configuration of P{kj, kj-^i), while the position of 
a' = a{kj^i+i') and (3' = f3{kj+i) in Fig. [H](a) determines the configuration 
of P{kj+i, kj+i+i'). As kj and kj^i run over all the values from —00 to 00, 
we find from Kronecker's theorem^ that every point in either Fig. [7|(b) or 
Fig. [H](a) is equally probable. However, the positions of (a,/3) and {a', (3') 
are completely correlated by the shift {a,b) in (15. lip , which is fixed as long 
as i and i' are unchanged. Thus, the joint probability for P{kj, /cj+i) to be in 
the mth configuration and P{kj+i, kj^i+i') to be in the m'th configuration 
is the area of the intersection of the two regions R{m) and R'{m'), and is 
denoted by 

Am,m'iiJ') = area(i?(m) n R\m')). (5.12) 

By superimposing Fig. [7|(b) on top of Fig. [H](a), we obtain Fig. [H](b). This 
figure gives the intersections of all the regions of Fig. [7|(b) with all the regions 
of Fig. [8|(a), and the joint probabilities can be read off as the areas of these 
intersections. 

For different values of i and i', we get different values of the width a 
given in (I5.1ip of the vertical slice in Fig. [8] (a) and also of the width b of the 
horizontal slice. After cutting and pasting the difference slices, the resulting 
figures are very different, so are the superimposed figures. Thus, the area of 
intersection Am,m'{^,^') depends on the choice of £ and £'. However, just as 
A{m) and A{m') do not depend on the shifts 7^, Am^m' also does not 
depend on these shifts. Moreover, it is easily seen that this joint probability 
Am,m' {£■, £') is not only the same for all the different regular pentagrids, but 
also the same for the different orientations of the parallelograms (i.e. different 
choices of j = 0, ■ ■ ■ , 4). 

We let i^t™-'"! denote the nth integer vector in the mth odd configuration 
for odd spins inside P, where = 1, ■ ■ ■ , N{m) and m = 1, ■ ■ ■ , 8. Similarly 
(jgnotes the ra'th integer vector of the m'th odd configuration of odd 
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spins inside P', with n' = 1, ■ ■ ■ , N{m') and m = 1, ■ ■ ■ , 8. From Tables 
Ill(a), Ill(b) and HVl the possible difference vectors ^i^I™'"! and ^i^l™''"'! 



for two spins in P and P' may be found. Adding, as in (14. 9 p and (15.81) . 
these to their corresponding reference integer vectors, with three of their 
components given in (14. 3p . (14. 5 p and (15. 6p . we obtain the two integer vectors 
K^iz) = for the spin in P and K°{z') = for the spin in P'. In 

(15. ip and (15. 3p we only need their difference 



K°iz') - K°{z) = i?!'"''"'] - Kt'"'"! = (4, 4, ■ ■ ■ , £4). 
We find from (gS]), ({52]) and ([53D that 



(5.13) 



■Jj+2 



[p-'e\ -£ + dK\ 



m,',n'] 
j+2 



7+3 



It is easy to see from (15. 7p and (15.110 that 

if ^ {a} < 1 



(5.14) 



5 



'3+2 



6j+4 ■- 



1 if 1-a ^ {a} < 1, 

if ^ {/3} < 1 - 6, 

1 if 1 - 6 ^ {/?} < 1, 



(5.15) 



which shows that in the four sectors of the unit square shown in Fig. [S](b), 
the 6i pairs are different. As a and b are functions of i and i' only, and the 



entries in Tables Ill(a), Ill(b) and IIVI are also independent of j, the results 



in (I5.14P are easily seen to be functions of i and and of {a} and {(3}, but 
they are independent of j. Therefore, we use the primed variables £,■ ■ ■ , £"" 
to denote these j-independent values of £j, ■ ■ ■ ,£j+4, i.e. 



9 pi 
'1 ^ 1 



(5.16) 



For different choices of j, [£o, ^1, ^2, ^3, ^4] is just a cyclic permutation of €. 

If we let the distance vector (£, £') between the two parallelograms P = 
P{kj,kj+i) and P' = P{kj + £, kj+i + £') in ([52]) be fixed, while letting both 
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kj and fcj+i vary from — oo to oo (which is equivalent to ^ oo), then the 
parallelograms P and P' each are in one of eight different odd configurations. 
Since the joint probability for P being in the mth configuration and P' in 
the m'th configuration is Am,m'{^,^') given in (15.121) . the double sum in (15. 2p 
can be rewritten as 

£,£' fcj,fcj+i K°{£)K°'{e') 
8 8 Nim)N(m') 

= xJ2J2Il AmM^^n E E u{K^-'-\K^"^''-'^). (5.17) 

e/' m=l m'=l n=l n'=l 

with Af = 5p, cf. Km . Using ([3ll]) and flSlTS]) . we find ([53]) becomes 

U(kM^ Kirn', n'l^ ^ COS [Re(g*ELo4C')] (^^')[.„A,-A]' (5-18) 

which is different for different j in view of (I5.16p . Since the correlation 
functions have the cyclic property shown in (13. 3p . and x°{q) can be evaluated 
by choosing parallelograms P and P' in (15.21) oriented in any one of the five 
directions (any choice of j), we can rewrite X°{qx,(ly) in a more symmetric 
way by expressing it as the sum over the five different orientations j, and 
then dividing the result by 5. This means, 

oo oo 

x°(q) = E E (5-19) 

^=— oo £'=—OD 

where (15.170 to (I5.16P and (13.30 are used to find 

^88 N{m) N{in') 

X°(q)«' = ^ E E ^m,m'{^j') E E ^(q'-^) (o"Cr')[£/'/",^"'/""], 
m=l m'=l n=l n'=l 

1 ^ 

c(q, ^) = 3 E cos Re [q*Cie + i'C + ^"C + ^'"C + ^""C')] • (5.20) 

It satisfies the following identities, 

r{^)-£,-£' = X°(q)£,^', r{^*)e,£ = X°(q)£,r, (5.21) 
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in which q* {(1x1—%) ■ The former identity in (15.211) is easily seen as a 
consequence of the reflection symmetry in the correlation function (crcr') = 
(cr'cr); the latter one is due to five-fold rotation and reflection symmetry^ 

In the actual calculation, because 5^ and 5^ generally differ in four sectors 
of the unit square, the contributions to the susceptibility from these different 
sectors are evaluated separately. 



5.2. Results 

To evaluate the wavevector-dependent susceptibility (15.191] , (I5.20p . we can 
compute the A^,™.' (^j the overlap area (15.121) and the €s from (15.141) . 
In Table [Til and Eq. (13. 2p in Section [3] we have expressed the pair-correlation 
function (aa')^^^^... ^^^^ in terms of Baxter's universal functions which can 
be evaluated using methods in our earlier work.''231'51ESl' 

Near the critical point k = 1 the leading asymptotic behavior of the pair- 
correlation function is the same Painleve III or V scaling function''^'^ as 
in the uniform rectangular lattice. Therefore, the scaling behavior of the 
central peak of in onr Penrose Ising model is also known*- '^^'^ to be 
the same as for the regular Ising model. 

More interesting is the incommensurate behavior of x(q) as a function 
of wavevector q and how it changes with temperature or k. At the critical 
point we expect to be a function that has everywhere dense 7/4-th 

power divergencies, but is locally integrable. It is nontrivial to show this 
directly, but this conclusion seems to impose itself as one approaches the 
critical point from either side. 

Since the correlation functions decay exponentially away from the critical 
point, we find that the are rapidly decreasing functions of £ and 

Putting terms of about the same order of magnitude together, we find 

oo I 

X°(q) = X°(q)o,o + $^5,, 5, = 2 ^ [x°(q).,n + X°(q*)£,n-i], (5.22) 
1=1 n=-e+i 

using both identities in (15.211) . We shall give density plots of several cases 
next, displaying the temperature dependence more clearly. However, it must 

^ In particular, one can start with the reflection symmetry about the direction of the A:j+3 
grid-line and its action on the parallelograms P{kj, kj^i) and P (fcj, fc ^-^j^). One arrives at 
{a'cr')[£ ii iri irir gii>i^ = {(T(j' ) [fj^i ^g"" J'" .i"] ■ Rcplaciug J 1 — j in (15 20|1 then completes the 
proof of the second identity in (15.211) . 



Q-Dependent Susceptibility on Penrose Tiles 



41 



be said that we can calculate x°(q) to very high precision in the cases shown, 
which fact is not clear from looking at these density plots. 

We shall give plots both above the critical temperature Tc, (/c> = k<l), 
and below (fc< = l/k <1). The value of modulus k corresponds to the 
row correlation length® 

^ = l/|arsinh(l/v^) - arsinh(v^)| (5.23) 

of the symmetric square-lattice Ising model for T > T^. For T < the true 
value of this row correlation length is ^/2 with ^ again given by fl5.23p . ^SMM 
At very low temperature, we only need to consider Se for very small i. For 
^ 2 the joint probabilities Am,m'{^,(^') can be easily evaluated by hand 
as it only involves the calculation of areas of triangles and rectangles. For 
fc< = .04847302, which corresponds to ^ ^ 1/2, we find Se < 10"^° for £ > 2. 
The density plot for l/x°(q) is shown in Fig. [9](a) for —An ^ qx,(ly ^ 47r 
where q = (g^, Qy)- We find ten-fold symmetry, corresponding to the five-fold 
symmetry of the Penrose tiling. 




-10 -5 5 10 -10 -5 5 10 



(a) (b) 

Fig. 9. Density plot of the q-dependent susceptibility showing l/x{<lx,qy) for the Z- 
invariant Ising model on a Pentagrid Lattice (Penrose Tiles) at very low temperature: (a) 
^ « 0.5, or A:< = .04847302; (b) ^ « 1, or fc< = .2363562, both for T < T^. 

For A;< = .2363562 (^ ^ 1), we find it necessary to consider all Si for 
£ ^ 4. As the temperature increases, larger and larger £'s are needed. To 
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evaluate the joint probability by hand is no longer feasible. To symbolically 
program the calculation for any values of £ and is highly nontrivial, as 
there are many different situations to take into account. It took us several 
months to sort out all cases, programming the calculation using Maple. 

A density plot for fc< = .2363562 is shown in Fig. [n](b). Plots of l/x°(q) 
for ^ ^ 4 and ^ ~ 8 are shown in Fig. [10] (a) and Fig. [in](b), together 
with corresponding plots for the dual cases with T > Tc in Fig. [TO](c) and 
Fig. [Tn](d). We can see clearly that the number of visible peaks increases as 
T ^ Tc and that this effect is more pronounced as Tc is approached from 
above. In Fig. [TTl a density plot for at ^ = 2 is given for — 167r ^ 

Qx, Qy ^ I6n, and we can already see some evidence for the quasiperiodic 
pattern of the q-dependent susceptibility in the full {q^, gy)-plane. 



5.3. The q-dependent susceptibility x'^Cq) 

Finally, we show that the q-dependent susceptibility x^{q) of the even sub- 



lattice is identical to of odd sublattice. From Tables Ill(a) and 



Ill(b) , we find there are 16 different even configurations for P{kj, fcj+i). The 
corresponding 16 regions are plotted in the unit square with {a} and {/3} 
along the axes in Fig. [T2](a). The regions for the 8 odd spin configurations 
of P{kj+1, kj+i+1) are also plotted with {a} and {[3} as axes in Fig. [T^(b). 
It is easy to see that after inverting one of the squares, the two figures are 
identical up to labeling. 

Furthermore, looking at Fig. [T2](b), we see that the five disjoint regions 
2, ■ ■ ■ ,6 become connected, if we impose periodic boundary conditions on 
the square. This "wrapping on a torus" is consistent with moving the slices 
as discussed below Eq. (15. lip . We next compare how the even regions in 
Fig. [12] (a) relate under the same periodic boundary conditions, to see if 
this is somehow true here also. As first examples we look at the two even 
configurations e(7) and e(13) and find that they have the same number of 
sites, and their difference vectors are related by 

= (0,-1,1) (5.24) 

For the other regions, we similarly find equal numbers of sites and 

6K^'^ = (0, -1, 1) + 6K'''\ = (0, -1, 1) + (5.25) 
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(c) (d) 

Fig. 10. Density plots of the q-dcpendcnt susceptibility, plotting l/x(q) versus and 
qy-. (a) fc< = .7018662 or ^ f« 4; (b) k< = .8379187 or ^ « 8, both for T < T^. Two 
corresponding plots for the dual models, with T > and identical values of fc> and ^, 
are given in (c) and (d). 
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Fig. 11. Density plot of the q-dependent susceptibility showing l/x(q) at A;< = .4912758 
or ^ w 2 and T < T^. 
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where m = 0, 1, 2 and n = 0, 1. 

The difference vectors of the odd configurations in Fig. [12] (b) are also 
related to those of the even configurations in Fig. [121(a), i.e. 

5K^^^ = (1, -1, 1) - 6K°'^, SK""'^ = (1, -1, 1) - 

= (0, 1, 0) - n = 0,---,5. (5.26) 

If the dependences of difference vectors on the values of {a} and {/?} are 
included in the equations, we may relate the even spins in P{kj, kj^i) with 
the odd spins in P{kj + 1, fcj+i + l) by 

6K^[{a}, {(3}] = (1, -1, 1) - (52, -62-6,, 6,) - 6K°[{-a}, {-/?}], (5.27) 

where 1 — {x} = {— x}, for x not an integer, and 

^2 = [1 - {«} + P''\ , ^4 = [1 - {P} + P''\ ■ (5.28) 

Consider the pentagrid with ■jj — > — 7^- and denote its parallelograms by 
P{kj, kj+i) such that 

a{-kj+i) = -a{k,+i), Pi-kj) = -Pikj). (5.29) 

Then it is easy to show that 

S2 + L-«(fc.+i)J = L«(l-%+i)J, <^4 + [-P{kj)\ = mi-k,)\. (5.30) 

Now, the integer vectors of even spins in P{kj, kj^i) given by (14.91) . (14. 3 P 
and (14. 5 p can be shown to relate to the integer vectors of the odd spins in 
Pil-k„l-k,+^) by 

i?^(e') = (l,l,l,l,l)-ie°(6). (5.31) 

This equation is consistent with the fact that the index of the odd spins 
is either 1 or 3, while the index of the even spins is either 2 or 4. Since 
the even spins in the original pentagrid are related to the odd spins in a 
different pentagrid, and the joint probabilities are independent of shifts, we 
have shown that the susceptibility of the even sublattice is identically the 
same as the one of the odd sublattice. 
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Fig. 12. (a) The sixteen even spin configurations of P{kj, fcj+i) plotted witli {a} and {/?} 
as the horizontal and vertical axis, (b) The eight odd spin configurations of P{kj+1, fcj+i+l) 
plotted with {a} and {/3}. 

6. Conclusions and Final Remarks 

In this paper we have presented a systematic way of evaluating the averaged 
pair-correlation function of a Z-invariant ferromagnetic Ising model with 
spins on half the sites of a Penrose tiling and Ising interactions across the 
diagonals of the rhombuses. 

Next, we have found that the q-dependent susceptibility of this model is a 
superposition of incommensurate everywhere- dense peaks, though not many 
peaks are visible at temperatures very far away from T^. For T < these 
peaks add a diffuse background to the Bragg peaks due to the spontaneous 
magnetization. Since the in (15.221) consists of Ai terms of the same order 
of magnitude, we compare their contributions. We find that the number 
of peaks of Se increases as £ increases, but that the numbers at fixed i are 
almost independent of temperature, even though the magnitudes of the peaks 
change as the temperature varies. 

As T — >■ Tc, the correlations decay more and more slowly, so the S^s 
increase, and more and more of these Sis are to be included in the numerical 
evaluation of the q-dependent susceptibility, which accounts for the ever- 
increasing number of peaks. This is unlike the behavior of the Fibonacci Ising 
models, considered earlier, ''^'^ where the ferromagnetic aperiodic Fibonacci 
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lattice behaves almost like the regular Ising model. 

Moreover, the q-dependent susceptibility is not a periodic function of 
or Qy. This behavior is different from that of aperiodic models defined on 
regular lattices.''^ This is because, when the lattice is aperiodic, we cannot 
separate the average of the correlation functions from the exponential (or 
cosine) terms, which contain the information about the lattice structure, as 
can be seen from (15.191) and (15.201) . 

At Tc, the q-dependent susceptibility has everywhere- dense divergences 
with the Ising exponent 7/4, but is still locally integrable. Away from Tc 
the x(q) is a continuous function. The Bragg peaks below Tc form a set of 
everywhere-dense Dirac delta functions of various strengths, but their sum 
is also locally integrable. These are strange objects and de Bruijn initiated 
their mathematical study. ''^ 

One of the main results of the current paper is that it provides a new 
method for doing calculations of probabilities on Penrose tilings. In Sec- 
tion m the calculation of the joint probability of the configurations of two 
parallelograms on the pentagrid is reduced to linear programming. 

Penrose tilings may be obtained by projecting certain subsets of the 
lattice into the plane.® The frequencies of the different types of vertices 
are then given as areas in the orthogonal spaces.® This method, which is 
known as the cut-and-project method, has been apphed and generalized by 
many 

authors. raiSaiSZllMllMiminiESlEaElESlEi The equival ence of the projec- 
tion method and a generalized grid method has been demonstrated. ''^'^'^ 
The positions of the Bragg peaks have been worked out,^'^^ '''^ ^ ''^-"^ to- 
gether with the values of probabilities of local configurations. The Penrose 
tilings can even be obtained from projections in a four-dimensional root lat- 
tice. 

It would be interesting to obtain results for joint probabilities similar to 
ours also by cut-and-project methods. This would generahze the windowing 
method of Baake and Grimm. We have not pursued this here, as our 
model is defined in terms of rapidity lines on the pentagrid. 

Another possible generalization of our work is to consider Penrose tilings 
that are periodic in either one or both directions. Finite approximants to the 
aperiodic Penrose tiling have been constructed through periodic pentagrids 
or projection methods. ''^'^'^'^ 

The mathematician Robinson has brought to our attention an exercise 
in the book by Griinbaum and Shephard''^ where a tiling with Penrose 
rhombuses can be cut into patches, and then converted into an aperiodic 
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set of 24 Wang's tiles. We have found that each patch in the Penrose tihng 
described in the exercise is in fact the image of a parallelogram under the 
mapping (12.51) . Thus these 24 configurations of the parallelogram can be 
easily converted into Wang's tiles. 

Finally, in our previous studies''^ we have examined the g-dependent 
susceptibility x(q) of some quasiperiodic Ising models on the square lattice 
defined in terms of Fibonacci sequences. It may be of interest to study 
models based on other sequences such as the aperiodic sequences studied by 
de Bruijn''^ or Tracy.®' We may ask what effect this has on the mixed 
interaction cases and also if it makes any difference for purely ferromagnetic 
models. 
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